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Abstract. In present paper an analytical approach is used to derive a new exact solution of a problem of one 
dimensional adiabatic flow of anon planner and cylindrical strong shock wave propagating in a mixture of 
gas and particle when density ahead of shock front is assume to vary as power of the distance from the source 
of explosion. A complete investigation is made for non - planer case and for cylindrical case an expression for 
total energy is obtained also. 


1. Introduction 

The investigation of shock wave phenomenon in a dusty gas is important in many 
engineering and astrophysical problems. Assuming solid particles to behave as a pseudo - 
fluid Pai° has derived appropriate basic equations for gas particle mixture. Using Lie Group 
method Jena and Sharma ‘have obtained entire class of self similar solution in dusty gases. 
Work done by authors Vishwakarma", Vishwakarma and Nath”, Vishwakarma and 
Vishwakarma” , Pandey and Shukla’ and Pandey and vaish* can also be cited in this 
context. Such type of study has numerous applications in underground explosions also. In 
this context work of Mcqueen? , Nagayama‘* and Lamb’ et-al can be cited. Sing and Husain 
1 has investigated analytical study of strong non planer shock in magnetogasdynamics. 
Ray’ has obtained an exact solution of a Spherical Blast wave under terrestrial conditions. 
Pandey and Singh® have studied non self similar solution of shock -wave in gas particle 
mixture. In present paper an analytical approach is used to derive a new exact solution of a 
problem of one dimensional adiabatic flow of anon planner and cylindrical strong shock 
wave propagating in a mixture of gas and particle when density ahead of shock front is 
assumed to vary as power of the distance from the source of explosion. A complete 
investigation is made for non-planer case and an expression for the total energy is 
obtained for spherical case also. 


2. Basic Equations 
Basic equation describing a planer, cylindrical and spherical motion for two phase flow of 
gas-particle mixture is given by Pai’ 
op A(pu) | ght 26, (2.1) 
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Ou r Ou 4 1 Op = 0, (2.2) 
Ot Or por 

gee P (2 ,uæ)-0, (2.3) 
Ot Or p°\ ôt Or 


where a= 0, 1 and 2 for planer, cylindrical, spherical symmetrical cases respectively and 
u, p,p and E are resultant velocity, pressure, density and internal energy per unit mass of 
the mixture. 

Here p=p,+p,, p,and p, being partial pressure of solid particles and gas respectively, 


p, =(1—-z)p and p,=zp, 
z being particle volume fraction. 
p=2p,, +(1-z)p,, 
Py and p, being species density of solid particles and partial density of gas, r is the single 
spatial co-ordinate being either axial in flows with planer geometry, or radial in 


cylindrically symmetric flows. 
The internal energy per unit mass is given by - 








p= lop (2.4) 
((-1)p 
h K, (2.5) 
wnere = 7 
i (I-K, )]G+K, 
G -Pr 
P; 


K „being mass concentration of pseudo-fluid of solid particles and is defined as 


Py. pa 
K, ===>. (2.6) 
ep P 
In equilibrium flow K,, is a constant in the whole flow field, therefore, from equation (2.6) 


~ = Constant. (2.7) 
p 
Thus equation of state becomes 
l-K 








= (egy deg) 
where R,, =(1—K,)R, and Ry may be considered as an effective gas constant of the 


mixture . 
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Here we have followed the analysis of two phase flows considered by Pai’ in which p 
is total pressure, z is the particle volume fraction and T is the ratio of specific heat defined 
by, 








= 2.9 
Cym 4 1+66 29) 
c K c 
VS ô=, ==> . (2.10) 
Cy l-K, Cy 
We can write 
c 





Cy 
Let 2 L, 
Cp 
then 
fri 
Y 
Thus equation (2.9) becomes 
T= yru) (2.11) 
(1+ uy) 


where c, and c, being specific heat at constant pressure and constant volume and 


C=C, ey; (2.12) 


p 


C,, being effective heat at constant volume for solid particle and c,is the specific heat of 


solid particles due to internal degree of freedom. 
Substituting value of E from equation (2.4) in equation (2.3) we get 


Ot Or p(U-z)\ ôt Or 





(2.13) 


Let, r = y (t) be the strong shock with the shock speed U = A propagating in to the 
t 
medium characterized by 
P=p(r),u=0, p=p,(r). (2.14) 
The conservation relations for mass, momentum and energy in present problem can be 
written as 
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Po = PU -u) | (2.15) 


Pot PU? =ptp(U-uy (2.16) 





Po l Pol- zo) ; U? _ Pı Pe) lig uy ; (2.17) 
Po p(T —1) 2 Pı oT -1) 2 

















Eee (2.18) 
Po P 
From equation (2.15) to (2.18), we get 
2 2 
ie 20120) U í a) (2.19) 
+ 
2(1-z,)eU?|, T-1143 
= 1 i 2.20 
: r+1 2r U? Coe 
P(E +1) 
pE Tr—1+2zo : (2.21) 
2 A 
1 2 1 
Terz ( U ) (1- zo) 
where 
(2.22) 


K > 
oe U-K, +K, 
andG, is the ratio of density of the solid particle to the initial density of the gas, p, and p, 
are respectively the pressure and density of the undisturbed atmosphere and A, is the 
T py l 
Py (l—Zo) 


Since the initial energy input E, of explosion is very large, the shock speed U >> A, 


sound velocity, given by Aj = 


A 
so that T — 0 in case of strong shock limit. 


Thus we have following strong shock conditions 


+1 
eee ee 2.23 
pP Oz, (2.23) 
„n-A -2)y, (2.24) 
+1 
2(1-zo) 2 
=———~ 9, U~. 2.25 
P Tel Po ( ) 


Assuming that p, varies as inverse power of radial distance from source of explosion, we 
have 


Po EPK] (2.26) 
where 
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0, are constants. 


The total energy in this case is given i - 


E, = (4K, "Te z) jfa =a r“dr. (2.27) 





3. Analytical solution 
With the help of equation (2.23) and equation (2.24) equation (2.25) becomes 








1| T-1+2z 2 
=—| ——— | pu’. 3.1 
f l- zo Je oe 
Differentiating (3.1) with respect to r and then putting value in equation (2.2) we get 
1} T-1+2 
u,, +uu,,.+ 2 | 2 Op + 2up 2] =0. (3.2) 
p| 2(1-z,) or Or 
With help of equations (2.18) , (2.23) , (2.4) and (3.1) equation (2.3) reduces into 
Pi 2 
u,, +UU,, + Mli pu +a% |=0. (3.3) 
2(1— zo) Or r 
Subtracting equation (3.2) from (3.3) and integrating it with respect to r, we get 
ur“ = é(t), (3.4) 


where ¢(t) is an arbitrary function of integration. 


Differentiating equation (3.4) with respect t and r then using equation (2.1) we get 


2 
Gig Ope OE aks, (3.5) 
r eb 
Subtracting equation (3.5) from equation (3.3) and after simplification , we have 


following linear differential equation 





i al Ee alow) Ee (3.5a) 
+1 r +1 Jé 
On solving this equation we get 
gee ee. (3.6) 
é Ot 
where 
2(1- zo) 





~ (P+ +(0+3)—2z,0 


5/9 


6/9 
KANTI PANDEY AND DEWKI NANDAN TEWARI 


0 
Differentiating equation (3.6) with respect to t and r and putting values of Spr we get 


2 2 
g A zi (3.6a) 
Ot é ôt 
This is a second order linear differential equation whose solution is given by 
G= $0 (3.7) 


where 


— (P+1)+(F+3)-2a 2 
— (P4i4+(P-lat+22,@ 
Using the jump condition (2.23) to (2.25) and equations (3.6) to (3.8) and equation (2.26) 
we have following relations 





T+ 


pero (3.8) 
(+1)+(U-l)a+2z,a 
5 
(T-3)a -(T +1)-4azo 
+1 
u=ndrt' =(1- z, Jort 
gorni (3.10) 





oO= 


Se (3.9) 





E, ory 


Using equation (2.27) and equation (3.10) the analytical expression for the total energy is 
given by 


22a Hae = Ds (3.11) 
+K, l-z\P-1) o Qa-2) 





4. Result and discussion 
In present article following Singh and Husian"’ , we have obtained a new exact solution of 
one dimensional adiabatic flow of non planer strong shock wave propagating in a mixture 
of gas and particle when density ahead of shock front is assumed to vary as power of the 
distance from the source of explosion. Variation of pressure , velocity and density is given 
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through figure (1,2,3) for z =.024,T =.97,@=1.97. Variation of energy is given by figure 
(4). 
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Fig2- variation of velocity for cylindrical and spherical symmetry. 
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Fig4- variation of energy for spherical symmetry. 
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